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Abstract. We continue the study of the effects of a strong magnetic field on the temperature distribution in 
the crust of a magnetized neutron star (NS) and its impact on the observable surface temperature. Extending 
the approach initiated in Geppert et al. we consider more complex and, hence, more realistic, magnetic 

field structures but still restrict ourselves to axisymmetric configurations. We put special emphasis on the heat 
blanketing effect of a toroidal field component. We show that asymmetric temperature distributions can occur and 
a crustal field consisting of dipolar poloidal and toroidal components will cause one polar spot to be larger than 
the opposing one. These two warm regions can be separated by an extended cold equatorial belt. As an example 
we present an internal magnetic field structure which can explain, assuming local blackbody emission, both the 
X-ray and optical spectra of the isolated NS RXJ 1856-3754, the hot polar regions dominating the X-ray flux and 
the equatorial belt contributing predominantly to the optical emission. We investigate the effects of the resulting 
surface temperature profiles on the observable lightcurve which an isolated thermally emitting NS would produce 
for different field geometries. The lightcurves will be both qualitatively (deviations from sinusoidal shape) and 
quantitatively (larger pulsed fraction for the same observational geometry) different from those of a NS with an 
isothermal crust. This opens the possibility to determine the internal magnetic field strengths and structures in 
NSs by modeling their X-ray lightcurves and spectra. The striking similarities of our model calculations with the 
observed spectra and pulse profiles of isolated thermally emitting NSs is an indication for the existence of strong 
magnetic field components maintained by crustal currents. 
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1. Introduction 

It has been understood since Greenstein & Hartke 
that in presence of a sufficiently strong magnetic field, 
> 10^° G, the surface temperature of a neutron star (NS) 
will not be uniform as is expected in the unmagnetized 
case. Anisotropy of heat transport, caused by both clas¬ 
sical and quantum magnetic field effects, in the thin low 
density (p < 10^° g cm“^) upper layer, the so called en- 
velope, results in a strongly reduced conductivity in the 
direction perpendicular to the field and an enhanced one 
along the field. As a result, the regions around the mag¬ 
netic poles are expected to be significantly warmer than 
the regions around the magnetic equator. Page and 

Page & Sarmiento dMlIi applied the Greenstein & Hartke 
formula to explain the lightcurves of the isolated ther¬ 
mally emitting NSs PSR 0833-45 (Vela), PSR 0656-bl4, 
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PSR 06304-178 (Geminga) and PSR 1055-52, consider¬ 
ing dipolar field configurations, and also the addition 
of a quadrupolar component, and including the General 
Relativistic curvature effects on photon trajectories. Much 
work has been dedicated to study the effects of the mag¬ 
netic field on the properties of the NS envelope and crust 
(for reviews, see Yakovlev & Kaminker 110041 Ventura & 
Potekhin I2()()ll Lai 1200 l|l . Due to the shallowness of the 
envelope, ^ 100 m, heat transport can be treated in the 
plane parallel approximation as a one dimensional prob¬ 
lem in which the heat flux is purely radial with a, locally 
uniform, magnetic field having some arbitrary orientation 
with respect to the radial direction. The most recent stud¬ 
ies (Potekhin & Yakovlev 120011 henceforth ‘PYOl, and 
Potekhin et al. I2003|l have included the best up to date 
transport coefficients and equation of state: they showed 
some deviations from the simple Greenstein & Hartke for¬ 
mula and finally give us a reasonably accurate description 
of a magnetized NS envelope. 
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In presence of a sufficiently strong magnetic field, 
~ 10^^ — 10^^ G, the anisotropy of heat transport, simply 
due to the classical effect of Larmor rotation of the elec¬ 
trons, extends to much higher densities and can even be 
present within the whole crust. Recently, we have shown 
(Geppert et al. 120041 subsequently Paper I) that in cases 
where the field geometry in the crust is such that the 
meridional component of the field dominates over its radial 
component in a large part of the crust, as is the extreme 
case of a magnetic field entirely confined to the crust, the 
non uniformity of the temperature, previously considered 
to be restricted to the envelope, may actually extend to 
the whole crust. This modified crustal temperature results 
in a different surface temperature distribution than what 
the simple Greenstein & Hartke model predicted. In con¬ 
tradistinction, in case the dipolar poloidal field has its cur¬ 
rents localized only in the core of the star, i.e., the field in 
the crust is considered as the one of a vacuum dipole, the 
crust is practically isothermal and the non uniformity of 
the surface temperature is then only due to field effects in 
the envefope. This result, that the geometry of the mag¬ 
netic field in the interior of the NS leaves an observable 
imprint at the surface, potentially allows us to study the 
internal structure of the magnetic field through modelling 
of the spectra and pulse profile of thermally emitting NSs. 

There exist growing observational evidence that the 
anisotropy of the heat transport in the envelope alone, 
assuming an otherwise isothermal crust, can not explain 
the surface temperature distributions of some observed 
NSs. The seven X-ray dim isolated NSs (XDINSs), dubbed 
’’The Magnificent Seven” (for reviews, see, e.g., Haberl 
120041 IMOdbl 1200 , ’tIi are nearby isolated NSs, all discovered 
in the X-ray band where they show a thermal spectrum. 
As a set they share several properties: 1) low X-ray ab¬ 
sorptions imply small distances, < 200 -500 pcs, 2) high 
blackbody temperatures deduced from the fit of their X- 
ray spectra indicate ages of the order of 10® yrs, possibly 
confirmed by measurements of proper motions which allow 
for tentative identihcation of their birth sites in clouds of 
the Gould’s belt, 3) no detection of radio emission, 4) no 
detection of a hard X-ray tail, typical of magnetospheric 
emission in active pulsars, 5) no association with a su¬ 
pernova remnant, 6) detected spin periods (in 5 cases) 
between 3 to 11 s, 7) estimates of their surface dipolar 
field strength Bq between 1 — 10 x 10^® G, either from a 
measurement of the period derivative P (1 case) or the 
interpretation of broad absorption lines (in 5 cases) in the 
X-ray spectrum as being proton cyclotron lines, 8) optical 
broad band photometric detections (5 cases) can be in¬ 
terpreted as being the Rayleigh-Jeans tail of a blackbody. 
However, these optical data are well above the Rayleigh- 
Jeans tail of the blackbody detected in the X-ray (“opti¬ 
cal excess”) and indicate the presence of an extended cold 
component of much larger area than the warm component 
observed in X-ray, the latter having an emitting radius 
(~ 3 — 5 km) much smaller that the usually assumed ra¬ 
dius of a NS(~ 10 — 15 km). Schwope et al. H2()()5|l tried 


to fit the lightcurve of RBS 1223 and concluded that only 
a surface temperature profile with relatively small, about 
4-5 kms across, hot polar regions may explain the observa¬ 
tions. Pons et al. (mni and Triimper et al. (mm arrived 
qualitatively at the same conclusion when they fitted the 
combined X-ray and optical spectrum of RX J1856.5-3754. 
In both cases, the smallness of the hot region is much be¬ 
low what can be reached by considering anisotropic heat 
transport limited to only a thin envelope. There are, how¬ 
ever, not many possibilities to produce non uniform sur¬ 
face temperature distributions. Principially, rotation leads 
to an oblateness which enhances the thickness of the inso- 
lating outer layer at the equatorial belt making that region 
cooler than the poles (e.g. Geppert & Wiebicke ll98I)|l . This 
effect is, however, even for the fastest rotating millisecond 
pulsars, too small and a fortiori neglibigle for those iso¬ 
lated thermally emitting NSs whose rotational period is 
in the order of seconds. Another way to heat the polar 
regions is the bombardment by ultrarelativistic charged 
particles as can be expected in radiopulsars (see Gil et 
ai.mm- Such heating is mostly confined to the polar 
caps whose angular radius can be roughly estimated as 
9pc ^ {2'KRjcPy^^, R being the star’s radius, c the speed 
of light and P the pulsar’s rotational period, giving a po¬ 
lar cap radius Cpc 0.46 P~^/'^ km (for P measured in 
sec.), which is much smaller than the estimated size, a few 
kms, of the warm regions in the X-ray dim isolated NSs 
whose periods are in the range of 3 to 11 s. 

The previous Paper I was devoted to demonstrate that 
a strong magnetic held can have signihcant effects on the 
crustal temperature distribution. We solved the station¬ 
ary heat transport equation in Newtonian approximation 
and considered very simple, purely dipolar poloidal crustal 
held structures. Our goal here is to consider more realis¬ 
tic models in which the currents maintaining the poloidal 
magnetic held are distributed between the core and the 
crust and we also consider the possible presence of a 
strong toroidal component in the crust. We will still con¬ 
sider only axisymmetric, dipolar, magnetic held conhgura- 
tion. Moreover, we now perform a wholly general relativis¬ 
tic formulation of the heat transport and energy balance 
equations. 

There exist also theoretical reasons that the presence 
of a strong toroidal held is very likely. In a series of papers 
Markay & Tayler 119731 (and references therein) showed 
that a stable magnetic held conhguration needs a coex¬ 
istence of poloidal and toroidal held components, having 
approximately the same strength. Recently, Braithwaite 
& Spruit (1200411 considered the stability of pure poloidal 
held geometries numerically in Ap stars and in magnetic 
white dwarfs and arrived at the same conclusion. 

Little is known about the magnetic held structure in 
NSs which is very likely determined by processes dur¬ 
ing the proto-NS phase and/or in a relatively short pe¬ 
riod after that epoch. A proto-NS dynamo (Thompson & 
Duncan ri993ll is unlikely to generate purely poloidal helds 
while diherential rotation will easily wrap any poloidal 
held and generate strong toroidal components (Kluzniak 
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& Ruderman llDDXI Wheeler, Meier, & Wilson E()()2II . The 
magneto-rotational instability (Balbus & Hawley I1991II 
also most certainly acts in proto-NSs (Akiyama et al. 

and results in toroidal fields from differential ro¬ 
tation (Balbus & Hawley 1199X1) . It is not yet explored 
whether the onset of supercondctivity in the core and 
the relatively fast crystallization of the largest part of the 
crust allows the field to relax into a stable (force free) 
state. However, it seems realistic to consider the effect 
of magnetic field configurations which consist of poloidal 
AND toroidal crustal components as well as of a star cen¬ 
tered poloidal one. The relative strengths as well as the 
positions of the crustal components will depend on the - 
certainly varying - creation processes of those fields and 
will therefore be varied in reasonable limits. 

The paper is organized as follows: In the next sec¬ 
tion, §12 we define the axisymmetric field structures and 
consider the resulting 2D heat transport and energy bal¬ 
ance equations. The next section, §0 presents our results 
for the crustal temperature distributions resulting from 
a large set of different field geometries and § 0 studies 
the observable effects of the resulting surface temperature 
distributions, assuming simple isotropic blackbody emis¬ 
sion. As an example we present a crustal magnetic field 
structure which allows us to reproduce the two tempera¬ 
ture spectral fits of RX JI856.5-3754. A discussion of our 
results and conclusions are the objects of § 0 


2. Heat transport and magnetic field 

We will assume that the heat transport through the core 
(pcore ~ 1-6 • lO^'^g cm“^) is practically unaffected by the 
magnetic field. This assumption is justified by the very 
likely occurence of proton superconductivity which, in the 
usually assumed case that it is a type H superconductor, 
will confine the magnetic flux to fluxoids of microscopic 
size leaving the bulk of the core matter devoid of magnetic 
field As a consequence we consider the stellar core to 
be isothermal and we only need to define the magnetic 
field structure and calculate the temperature distribution 
in the crust. Moreover, we will still, in the present work, 
restrict ourselves to axisymmetric configurations. 

2.1. St ret u re of the Magnetic Field 

The assumption of axial symmetry suggests a decomposi¬ 
tion of the NS magnetic field in a poloidal and a toroidal 
component (we refer the reader to Radler l^OOOI for a gen¬ 
eral presentation of the decomposition of magnetic fields 
into poloidal and toroidal components) and we will restrict 
ourselves to the lowest order terms, i.e., dipolar fields. The 
boundary conditions at the NS surface are such that the 
internal field matches smoothly an external potential field; 

^ In case the protons form a type I superconductor, as 
has been proposed by Link (1200311 and Buckley, Metliski & 
Zhitnisky 1200411 . this assumption would have to be revised. 


an observer above the NS surface would see only a dipolar 
poloidal field, whatever internal field structure is present. 

Basically we will consider three components of the NS 
magnetic field: 


1. The core field is a dipolar field having no cur¬ 

rents in the crust and is maintained by axisymmetric cur¬ 
rents in the core. Then a vacuum-like dipolar field will 
penetrate the crust with components, in spherical coordi- 
inates (r, 9, ip) oriented along the symmetry axis. 


Dcore Dcore 
— ^0 


cos 9 


( 1 ) 


DCore 

^9 


Hg”® sin6» 


( 2 ) 


= 0 , (3) 

where H™"'® denotes the value of this component of the 
magnetic field at the magnetic pole on the surface of the 
star and x = r/R^s the dimensionless radial coordinate, 
normalized on the NS radius Rns- This description of the 
core field in the crust assumes that it feels the crustal shell 
above the crust-core boundary as vacuum. 

2. The crustal poloidal magnetic field is main¬ 

tained by toroidal currents circulating in the crust. As for 
the core field we will consider here only its largest scale 
component, which continues beyond the NS surface as a 
dipolar vacuum field and can be conveniently described 
in terms of the (possibly time dependent) Stokes stream 
function S = S(r,t). 

The field components are then expressed as 


TDCrUSt 


2 cos 9 


S{r, t) = 


cos 9 


s{x, f) 


( 4 ) 


Dcmst sin 6>95'(r,f) _ sm9 ds{xR) 

® r dr 2 X dx 

H®™®* = 0 (6) 

where is - as for B™’’® - the strength of this field 

component at the magnetic pole. We have introduced the 
normalized function s = 2S/Bo/Rlfg, which must satisfy 
the boundary conditions: s = 0 at the crust-core boundary 
and s = 1 with ds/dx = — 1 at the surface. 

3. The crustal toroidal magnetic field B*®’' is main¬ 
tained by poloidal currents circulating in the crust. In an 
axisymmetric configuration it vanishes at the NS surface 
and we will assume it does not penetrate the core. 

In these condition has vanishing components in 
the r and 9 directions and, applying the same formalism 
as led to the representation of the poloidal field by the 
scalar quantity s (see e.g.: Radler & Geppert ITTMI) , we 
can express the dipolar component of the toroidal field 
via the scalar t(r, t) by 


B*°r = —T(r, t) = B*®® — 


i(x, t) 


( 7 ) 


r 


X 
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and t must vanish at the surface and at the crust-core 
boundary while its maximum value within the crust is 
fixed to 1. Notice that this component is symmetric with 
respect to the equatorial plane while the dipolar poloidal 
fields are antisymmetric and, hence, the total magnetic 
field is asymmetric. To maintain an antisymmetric field 
configuration, to a dipolar poloidal field one would have 
to add a quadrupolar toroidal cx sin 0 cos 0, with all 
toroidal components with odd multipolarity vanishing. 

Determination of the spin-down of a NS can give us an 
estimate of the strength of the dipolar poloidal field at the 
magnetic pole, Bq. In our models Bq = Bg™®*-l-iJQ”®: this 
gives us no information about either the possible values of 
and the relative strengths of B™"'® and or 

about the shape of the two functions s{x) and t{x). There 
is thus a large degree of arbitrariness in the choice of the 
field structure, even within our restriction to the dipolar 
terms, and we will consequently consider many differente 
cases. We choose for s(r) one of the models of Page et al. 
I12()()()ll at an age t ~ 10® yrs: in these time evolution mod¬ 
els of crustal field, the field was not allowed to penetrate 
the stellar core due to proton superconductivity and the 
crustal currents sustaining the field were slowly migrating 
toward the high conductivity layers near the crust-core 
boundary. This gives us some confidence that our choice 
of s(a;) is not too unrealistic^. For the function t{x) we use 
a Gaussian-like shape with a maximum at 


= Xc 


aAxc 


( 8 ) 


where Xcore and Aa;crust are the normalized core radius and 
crust thickness, resp., and consider three values of a: 0.2, 
0.5, and 0.7. The normalized Stokes function s(a:) and the 
three different t(x) functions we will consider are plotted 
in Fig. 1 and the crustal field lines of the fields are pre¬ 
sented in the panels of Fig. 2. It is important to notice that 
the maximum values of is ~ I.3B™"'®, 

reached at the crust-core boundary, the maximum value of 
Btor jg Bq°’^, in the crust at the point x = x^ax, while the 
maximum value of is about 13 times from 

its 0 component (Eq. EJ, due to the fact that s{x) must 
decrease from 1 at the surface to 0 at the crust-core bound¬ 
ary and hence have a large derivative ds/dx ^ As/ AcCcrust 
with Aa:crust ^ 0.1, as shown in Fig. 1 


2.2. The two-dimensional heat transport 

The thermal evolution of the crust is determined by the 
energy balance equation which has, in axial symmetry and 
taking into account general relativistic effects, the follow- 


^ Notice, however, that the models of Page et al. used 

a ID numerical code and hence no 8 dependance of the tem¬ 
perature in the crust, which is in contradiction with the results 
of §0 


ing form: 
e 


.-A Q 


dr 


o2$ 


(r^ Fr + 


d 


r sin 9 dO 


{sinOFe) = 


at 


( 9 ) 


where T is the local temperature, e*, are the redshift 
and length correction factor, Fr and Fq are the local radial 
and meridional components of the heat flux and r and 9 
the local coordinates. Qi, and are the neutrino emis- 
sivity and specific heat, respectively, per unit volume. We 
will study stationary configurations and neglect neutrino 
energy losses, i.e., replace the right-hand side of Eq. El by 
zero, and rewrite this equation as 


dx 


{X^ Fr) 


1 


d 


X sin 9 d9 


[sin 9 Fg) = 0 


( 10 ) 


where d/dx = e~^djdx and iA.e = e^'^Frfi/RNS 

In regions below the envelope, i.e. in the crust, the 
quantized motion of electrons transverse to the magnetic 
field lines doesn’t play any role for the magnetic modifica¬ 
tion of the heat transport. The magnetic field acts via the 
Larmor rotation of the electrons, the relevant heat car¬ 
riers, impeding on the heat transfer perpendicular to the 
crustal field. Introducing, as in Paper I, the components of 
the temperature gradient parallel and perpendicular to the 
magnetic field as well as the corresponding components of 
the heat conductivity tensor we arrive at the heat flux 
vector 

[V(e^T) + b (V(e*T) • 6) + 6 x V(e^r)] (11) 


whose radial and meridional components are 


Fr = -Yi - - 



d9 2x'^ dx 


X3^-^sm6' 

x"' d9 


9^-0 a 2^ 

Sin 9 cos 9 —— —7 cos 9 


dx 


( 12 ) 


Fg = - 


1 dT 

xi--^ 
X d9 


.. dT s ds . ^ r, dT 1 fds\ . 2 ^ 


~tdT. 

- Y3-^sm6'. 

X dx 

where we have written T = e®T and Xi = e^Xi- 
The heat conductivity coefficients are: 


(13) 


Xi = 


Ko 


1 -I- (WbT)^ 


X2 = 


^0 (^Bgr)^ 

1 -I- (wbt)^ 


X3 = 


Ko(wB‘r) 

1 + (WbT)2^’ 

(14) 
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Note that the Larmor-frequencies ujgp,t in Eq.^^are calcu¬ 
lated with Bl = -b in X 2 and with 

for the “toroidal” heat conductivity coefficient X 3 . The 
Larmor-frequency in the denominator of these coefficients 
is calculated by using the exact field strength at each 
point: ^2 = Of + 9)^ + 0)^. 

For a given magnetic field structure equation IHl is 
solved with the heat flux components of equations El and 
1 1 31 numerically using the same code which has been briefly 
described in Paper I until a stationary solution is found. 
The neutrino emissivity in the crust can be neglected for 
isolated NSs of that age. The inner boundary condition, at 
the crust-core boundary is T(a; = Xcore, 9) = Tcore, be. the 
core is assumed to be isothermal. For the outer boundary 
condition the argumentation given in section 2.4. of Paper 
I remains valid and we proceed in the same way: we stop 
our integration at a density pb = 10 ^° g cm“^ and, at the 
corresponding a;b and for each 0 , the radial flux is given 
by Fr{xh,9) = (JsbTs where Tg is the surface (effective) 
temperature corresponding to the boundary temperature 
T{xh,d) obtained with the “Tb — Tg” resulting from the 
magnetized envelope models of PYOl. 

We use the same input physics as in Paper I for calcu¬ 
lating the electron relaxation time t and the same chem¬ 
ical composition of the crust. For the core of the star we 
here used the recent models of Akmal, Pandharipande & 
Ravenhall dMHI. We consider only a 1.4 Mq NS, with a 
radius Rns = 11-4 km, R°^ = Rns ' = 14.2 km, and 

a^crust = 0.914, Axcrust = 0.086. 

The essence of the effect of a strong magnetic field is 
that the heat flux F is forced to be almost aligned with 
the local field S when (wbt)^ 1 since then the com¬ 
ponent of the thermal conductivity tensor k parallel to 
B is K|| = kq while in the perpendicular directions it is 
K± = Ko/(l + (ojbt)^) ^ K||- Values of WbT were plotted 
in Fig 2 of Paper 1. Notice that we do not calculate By 
axial symmetry this component is independent of <p but 
certainly not equal to zero, in spite of having dT / d(p = 0^. 
Since, for strong fields, heat essentially flows along the field 
lines, when is dominant will also be much larger 
than Fg and Fr and produce a winding of the heat flow 
around the symmetry axis: F follows the shortest possi¬ 
ble paths with the highest possible conductivity and this 
winding effectively acts as a heat blanket. 

3. Crustal temperature profiles 

In this section we will examine the temperature distri¬ 
butions resulting from a set of representative magnetic 
field configurations in which the three dipolar components 
^crust^ Fc°"'e, and have different relative stengths. 
We are faced with many parameters since beside the 
strengths of each component we have also the essentially 
unknown locations of the crustal currents sustaining the 

® This is a strong limitation of onr simplifying assnmption of 
axial symmetry and can only be removed by using a 3D code 
which, hopefully, will be studied in the near future. 


crustal fields, i.e., the two functions s and t of Fq. 0 H 
and [71 The cases of purely poloidal fields, i.e., B*” = 0, 
with only or but not both components su¬ 

perposed, were studied in Paper I and we will describe 
here a large set of models which, we hope, will allow us 
to get a good understanding of the joint effects of these, 
still simplified axisymmetric dipolar, multicomponent field 
configurations. 

We consider first the case where the whole magnetic 
flux is produced in the crust and the field does not pene¬ 
trate the core, i.e., B = _|_ ^tor ^core = 

3 and 4 show the crustal temperature distributions, and 
the resulting surface temperatures, for two different loca¬ 
tions of the toroidal components, with Xma.x = 0.7 and 0.2, 
respectively, and in each case for three different strengths 
of Panels a) and b) of Fig. 3 show that, given the 

high Tcore, a Fg™®* of 10^^ or 10^^ G is not strong enough 
to affect the heat flow, in agreement with our previous re¬ 
sults of Paper I: the inner part of the crust is essentially 
isothermal and it is the strong which act as a heat 
barrier, through the winding of the field lines in the (f di¬ 
rection. With a very strong however, as in panel c) 

of Fig. 3, the thermal insulation is then practically entirely 
produced by this component and is located too far 
out to play any noticable role. In Fig. 4 we illustrate the 
same cases of but with F*®® now located at higher 

densities: for Fg™®* = 10*^ and 10*^ G, panels a) and b), 
we see again that the thermal insulation is provided by the 
toroidal component but with the new effect that heat can 
flow “around” the toroidal field, from the magnetic axis 
toward the equator, in the outer part of the crust and re¬ 
sult in much smaller meridional temperature gradients in 
this region. (This flow of heat “around” the toroidal field 
can also be seen in Fig. 3 but is much less efficient.) The 
case of a very strong F®''®‘®*, panel c) of Fig. 4, with this 
inner toroidal component is almost identical to the similar 
case of panel c) of Fig. 3. 

In the next two Figures 5 and 6 we explore the effect 
of the distribution of the poloidal flux between the crust 
and the core, together with a toroidal component, i.e., 
we consider the three components F®™®*, F®®’™, and F*°® 
together. We restrict ourselves to a F*®® whose maximum 
is located deep inside the crust, i.e., a = 0.2 or 0.5 in 
Fq.m which, as seen above, is not as effective as when 
located further out but is possibly more realistic since one 
expects currents to migrate toward the highest electrical 
conductivity layers (Page et al. I2000|l . In panels a) and 
b) of Fig. 5 the crust and the core contribute equally to 
the poloidal field and we see that when it is much weaker 
than the toroidal component, panel a), the latter acts as 
a heat barrier while in panel b) the poloidal component is 
sufficiently strong to allow heat propagation in the radial 
direction, superposed of course to the winding flow in the 
(p direction, along the field lines and only the region within 
about 30° above and below the equator where the field 
lines are not connected to the core are slightly colder. In 
case the poloidal flux comes almost entirely from the core, 
Fig. 5 panel d), the insulating effect of the toroidal field 
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almost completely disappear. Assuming equipartition of 
the poloidal flux between the core and crust components, 
i.e. ^ 10% coming from the crustal comonent, as in Fig. 
5 panel c), the result is very similar to panel c) despite of 
the fact that in this case the maximum of is about 

1.3 X 10^^ G (from its 9 component, see Fig. 1), i.e., larger 
than B™"®. 

We consider now an even stronger B^°'^ in Fig. 6 and 
see that its insulating effect is again very strong, due to 
the presence of the (wbt)^ in the denominators of Fn. 1141 
which is now increased by almost one order of magnitude. 
The results are perfectly in line with our description of 
the other previous cases: 1) the two models where the 
poloidal flux is crust dominated exhibit larger insulating 
effects than the corresponging models where the flux is 
core dominated, see panels a) vs. c) and panels b) vs. d), 
and 2) models where the toroidal field is located deeper 
inside the crust are less insulating than when located fur¬ 
ther out, see panels b) vs. a) and panels d) vs. c), since 
in the former cases heat can more easily flow arround the 
toroidal component back from the symmetry axis toward 
the equator. 

A noticeable general feature of all our results is the 
asymmetry between the two magnetic hemispheres, re¬ 
sulting from the asymmetry of the total field B, since 
^crust ^core gj.g anti-symmetric while 5*°“' is sym¬ 
metric. In all cases where a non-isothermality is apparent 
the warm region close to the symmetry axis is broader 
in the lower hemisphere than in the upper one. In cases 
where the poloidal component is almost comparable to the 
toroidal one, as in panels c)’s of Fig. 3 and 4, the asym¬ 
metry is barely distinguishible but in all other cases it is 
clearly visible. 

The lower panels of the four Figures 3 to 6 show the 
resulting surface temperature profiles T(6*), 6 being the co¬ 
latitude measured from the symmetry axis, scaled to the 
temperature at the magnetic pole 0 = 0, and compared, 
as dotted lines, to the naive cos^/^(0b) scaling (which it¬ 
self is a simplification of the Greenstein & Hartke, 119831 
result). These allow to easily produce the temperature dis¬ 
tribution over the whole NS surface and generate compos¬ 
ite spectra and pulse prohles which we describe in the next 
section. 

4. Observational consequences for thermally 

emitting neutron stars 

The very distinct surface temperature distributions result¬ 
ing from non-isothermal crusts as calculated in the previ¬ 
ous section, and displayed in the lower panels of Figures 
3 to 6, have several immediate obervational consequences 
which we explore in this section. In all our simulations of 
spectra and pulse profiles we use simple isotropic black- 
body emission at each point of the stellar surface, with the 
local temperature T(9, ip), and include general relativistic 
light bending as described in Page dMB. 

In presence of a strong toroidal field in the crust, the 
channeling of heat toward the two polar regions results 


in the appearance of two hot spots of very reduced size. 
The Figure 7 shows five examples of surface temperature 
distributions and the resulting observable pulse profiles. 
Naturally, models with the smallest hot spots result in 
the highest pulsed fractions, Pf, with values above 30%, in 
contradistinction to the case of an isothermal crust which 
results in P/ ~ 5 %. 

The composite blackbody spectra resulting from the 
same five cases of Figure 7 are shown in Figure ??. 
Assumed distance to the stars have been adjusted to give 
the same maximum flux in the X-ray band, and thus very 
similar X-ray spectra: given this adjustment the difference 
between the relative areas of the hot and cold regions in 
the various cases result in differences in the predicted op¬ 
tical fluxes. Gomparison, by eye, of the surface tempera¬ 
ture plots, left panels of Figure 7, with the relative opti¬ 
cal fluxes shows a direct correlation between the relative 
size of the cold region with the optical flux. This is con¬ 
firmed by Figure 9 in which the percentage of the total 
surface area whose temperature is below a certain value is 
shown and exhibit the same ordering as the optical fluxes 
of Figure 8. This correlation can easily be made quantita¬ 
tive by noticing that most of the flux comes into X-ray and 
the latter is thus given by Fx — a{R/D)‘^asBT^ff, D being 
the star’s distance and a a correction factor for the inter¬ 
stellar absorption, where the effective temperature Tea is 
naturally defined by 


TeS = 




1/4 


(15) 


The optical flux, however, corresponding to the Rayleigh- 
Jeans tail of the spectrum is given by Fq — {R/D)‘^ATa-ve, 
A being a constant dependent on the energy range in¬ 
cluded in the ’’optical”, and the average temperature Tave 
being also naturally dehned as 


T = 


dn 

Att 


Ts{0,ip) 


One hence has the simple relationship 


Fo — Fx 


Ta, 


A 


rp4 


a (JSB 


(16) 


(17) 


The spectra plotted in Figure 8 considered different dis¬ 
tances D, adjusted to produce almost identical X-ray flux 
but with the same Nh and hence approximately the same 
a’s, and Eq. El allows to correctly predict the relative 
Fo’s, within 10%, using the values of (Tave,7eff) given in 
Figure 7. 

As a concrete example, we show in Figure 10 an at¬ 
tempt to model the spectrum of the isolated NS RX 
J1856.5-3754. Its X-ray spectrum is well fitted by a hot 
blackbody, with T°° ^ 6 — 7 x 10® K and an effective 
radius R^ ^ 4.5 km, while the optical data can be in¬ 
terpreted as the Rayleigh-Jeans tail of a cold blackbody 
component with T°° < 3.5 x 10® K and an effective radius 
>17 km, assuming a distance U = 117 pcs (Pons et 
al. 120021 Triimper et al. I2004|l which are plotted in this 
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figures as dashed lines. Using the field structure of panel 
c) of Figure 6 and the resulting surface temperature pro¬ 
file we obtain a reasonable ht, shown as a continuous line, 
to this two blackbodies description. 

5. Discussion and conclusions 

We have studied the impact of different crustal magnetic 
held conhgurations on the transport of heat within the 
crust of a NS and the resulting surface temperature distri¬ 
butions with their observable consequences. In a previous 
paper (Paper I) we had hrst considered the simplest cases 
of a poloidal dipolar held either entirely conhned to the 
crust (“crustal held”) or exclusively produced by currents 
in the core and hence having in the crust the same struc¬ 
ture as a vacuum dipole (“core held”). In the present work 
we have considered the effect of an additional (dipolar) 
toroidal component, located within the crust, with vari¬ 
ous combinations of core and crustal helds, and all three 
components having the same symmetry axis. 

Since we consider here only poloidal and toroidal dipo¬ 
lar components, due to the axial symmetry the toroidal 
held component is conhned to the interior of the NS and 
the external held geometry is always a purely poloidal 
dipolar one uniquely characterized by its strength at the 
magnetic pole, Bq. Therefore, the effects of various mag¬ 
netic held structures with the same Bq in the magneto¬ 
sphere, as e.g. pulsar properties and rotational evolution 
of the NS, are indistinguishable. However, these differ¬ 
ent inner geometries of the magnetic held result in very 
different surface temperature distributions which are po¬ 
tentially distinguishable through observations of thermal 
emission from isolated cooling NSs. 

We had shown in Paper I that a strong magnetic held 
practically channels the heat how along its held lines. We 
had found that a core held, with held lines essentially ra¬ 
dial in most of the crust, results in an almost isothermal 
crust and the surface temperature is hence only controlled 
by the anisotropic transport in the uppermost layers of the 
envelope (Greenstein & Hartke lT983l Page 1^23) exhibit¬ 
ing two symmetric extended warm regions, one in each 
hemisphere, and a cold belt along the magnetic equator. 
In contradistinction, a crustal held, having held lines with 
large meridional components, can produce strongly non- 
isothermal crusts and the heat how is much more con¬ 
centrated around the magnetic poles, resulting in a much 
wider cold equatorial belt. Our main result here is that the 
addition of a toroidal component, inside the crust, is able 
to add a very efhcient heat blanket which may force the 
heat to how within a narrow region along the polar axis 
and result in a cold region covering most of the stellar sur¬ 
face with two warm spots arround the two magnetic poles. 
However, we hnd that for this effect to occur it is essential 
to have a signihcant part of the poloidal held produced 
in the crust: the toroidal held in itself is not sufficient in 
case the held lines of the poloidal component are almost 
radial. 


The presence of two small warm regions separated by 
an extended cold belt has two inmediate observational 
consequences. The hrst is that the observable pulsed frac¬ 
tion in the X-ray band can be very large, above 30% 
assuming isotropic blackbody emission, and, second, the 
cold region’s emission contributes little to the X-ray hux 
but dominates the detectable hux in the optical range, ap¬ 
pearing as an “optical excess”. Notice, however, that addi¬ 
tion to the dipolar poloidal held of a quadrupolar compo¬ 
nent also allows, potentially, to reach high pulsed fractions 
(Page & Sarmiento llDDtill within the same isotropic black- 
body emission scheme and without any todoidal compo¬ 
nent. In this latter case, high Pf’s are reached when the 
quadrupole sufficiently deforms the held to push the two 
magnetic poles close to each other but without signih- 
cantly altering the star’s ehective and average temper¬ 
ature (see, e.g.. Figure 1 of Page & Sarmiento I199till . 
and hence producing no signihcant “optical excess”. The 
“Magnihcent Seven” briehy presented in the Introduction 
seem to be very good candidate NSs for which our re¬ 
sults may be relevent. Their estimated dipolar surface held 
strengths Bq are above 10^^ G, opening the possibility of 
strong anisotropic heat transport in their crust, and the 
“optical excess” observed in 5 cases is then a natural con¬ 
sequence of the large Hq’s in case their crusts harbor very 
strong poloidal helds. As we has seen, channeling of the 
heat hux toward two small polar regions requieres more¬ 
over that part of the poloidal hux be conhned to the crust: 
this can be achieved either by the action of some dynamo 
process which produced such a poloidal held component 
or/and by the expulsion of hux from the core through the 
action of the outward migration of the neutron superhuid 
huxoid resulting from the NS spin-down and the consecu¬ 
tive expulsion of the magnetic huxoid resulting from the 
core proton superconductor. This latter process is likely to 
have occured in the “Magnihcent Seven” given their un¬ 
usually long period, for their relatively young ages, which 
are very likely the results of fast spin-down due to their 
strong dipolar poloidal helds. However, our axisymmet- 
ric held conhgurations produce symmetric, but not sinu¬ 
soidal, light-curves and the superposition of a quadrupo¬ 
lar component may be necessary to produce precise hts of 
the observed pulse prohles of the seven XDINS (Zane & 
Turolla l!l()()5ll 

Results very similar to the ones presented here have 
been recently obtained by Perez-Azorin, Miralles & Pons 
who also conclude that a toroidal component local¬ 
ized within the crust lead to surface temperature distri¬ 
butions with very localized warm regions and an extended 
cold equatorial belt. These authors also took into account 
the transport of heat by phonons, which is unaffected 
by the magnetic held and thus reduces the anisotropic- 
ity, and the small quantum corrections present at high 
densities for very strong helds, two effects we did not con¬ 
sider. Moreover, they performed their heat transport cal¬ 
culations to lower densities than we have considered here 
and assumed that the upper layer is cut at a hnite den¬ 
sity due to solidihcation of the material while we simply 







Geppert, Kiiker & Page: Temperature distribution in magnetized neutron star crusts II 


assumed that the envelope extends into an atmosphere 
and “glued” our models with pre-existent envelope mod¬ 
els. Nevertheless, despite of these differences, the similar¬ 
ity of the results prove that the effect of a toroidal held is 
independent of the details and make our common results 
very robust s. 

Claiming that the observed optical-X-ray properties 
of the “Magnihcent Seven” prove the existence of strong 
toroidal helds in their crust may still be premature. Given 
our poor understanding of the emissive properties of the 
NS surface, being it either an atmosphere, a liquid or a 
solid depending on its chemical composition, temperature 
and magnetic held strength, one cannot exclude that the 
“optical excess” is achieved locally, i.e., is an intrinsic fea¬ 
ture of the spectrum produced at each point of the sur¬ 
face, instead of being a global effect due to strong sur¬ 
face temperature differences as we have obtained. In con¬ 
tradistinction, three very active pulsars exhibit very diher- 
ent properties (Vela: Romani, Kargaltsev & Pavlov 2005, 
120051 PSR B0656-I-14: Shibanov et al. 120051 and Geminga: 
Kargaltsev et al. 1200511 : a thermal spectrum is clearly seen 
in the X-ray but near & far UV observations give upper 
limits on its Rayleigh-Jeans tail which are in agreement or 
even below the extrapolations of the X-ray spectrum. The 
reasons for this diherent behavior is presently unknown 
and may be due to their weaker magnetic held, different 
surface chemical composition, and/or reprocessing of the 
thermal photons in the magnetosphere. 

Finally, an issue raised by our models is the time evo¬ 
lution and stability of the held conhgurations we have 
considered. We intend to perform coupled time-dependent 
cooling and held evolution calculations to try to eliminate 
at least the most unlikely conhgurations, i.e., either highly 
unstable conhgurations and/or conhgurations based on 
current distribution which will evolve fast and reajust or 
decay. Nevertheless we can speculate that a very strong 
crustal toroidal held could produce such strong tension 
that the crust should readjust itself, as is seen in magne- 
tars, and this may be the cause of the observed change 
in the spectrum and pulse prohle of RX J0720.4-3125 (de 
Vries et al. l2004ll . 
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6. Figure captions 

Fig. 1 : Gontinuous lines: normalized Stokes function s(a:) 
and its derivative (scaled by a factor 10) used in this work, 
Eq.4 and 5. Discontinuous lines: the three different nor¬ 
malized functions t{x) we consider for the toroidal held, 
Eq.7. 

Fig. 2: Qualitative sketch of the held geometries as con¬ 
sidered here. Left panel: held lines of the poloidal crustal 
held having outside the same dipolar potential held as a 
star centered held (right panel). Middle panel: isolines of a 


toroidal held with its maximum in the middle of the crust 
(a = 0.5). 

Fig. 3 : Thermal structure of the crust of a neutron star. 
(The radial scale of the crust is stretched by a factor of 
5 for clarity.) The magnetic held includes the and 

B^°'^ components and its structure is illustrated in panel d 
(held lines of the poloidal component and color coded iso¬ 
lines of the toroidal component normalized on The 

maximum of 5*°"' is located at Xmax = a^core 4-0.7 AcCcrust 
(see Eq. 8). Panels a,b, and c show the temperature dis¬ 
tribution for Rj” = IQi'^G and = lO^^, 10^3, and 

lO^^G, respectively. The color coded temperatures are nor¬ 
malized on the core temperature, Tcore = 6 • lO^K. The 
lower panels show the resulting surface temperature pro- 
hles T{9), once an extra insulating envelope has been glued 
to the core temperature distributions of the upper panels, 
and the dotted lines illustrate the same prohle when an 
isothermal crust is assumed. 

Fig. 4 : The same as in Fig. 3 but assuming a location of 
the maximum of the toroidal component much closer to 
the crust-core interface: Xmax = Xcore 4- 0.2 Axcrust (see 
Eq. 8) 

Fig. 5 : The same as in Fig. 4 but now including the three 
components and S*”. The magnetic held 

scales are in panel a: = IO^^q^ 5 ™re ^ 19120 , 

panel b: Rg™®* = lO^^G, Rg”® = lO^^G, panel c: R§™®* = 
IQi^G, Rg”® = lO^^G, and panel d: R§™®* = lO^iQ, 
= 1013 g. In all four cases Rj®’' = lO^^ G and 
^max — ^core 4- 0.2 Axcrust (see Eq. 8). 

Fig. 6: The same as Fig. 5 but assuming a stronger 
toroidal held with Rq®*' = 3 - lO^^G. In panels a and b the 
crustal poloidal held dominates the core held (Rq'’®®* = 
7.5 • IQi^G, Rg®"® = 2.5 • IO^^g) while in panels c and d 
the core held is dominant (Rg™®* = 2.5 • IO^^g, 5 ™re 
7.5T0^^G). In panels a and c the toroidal held has its max¬ 
imum in the middle of the crust, Xmax = Xcore 4-0.5 Axcrust 
(see Eq. 8), while in panels b and d it is located closer to 
the crust-core boundary Xmax = Xcore 4-0.2 Axcrust- 
Fig. 7 : Surface temperature distributions (left panels) in 
an area preserving representation with a color scale follow¬ 
ing the emited hux (oc T'^). Panels a to d use the internal 
held structures of the corresponding panels in Fig. 6 while 
panel e assumes an isothermal crust. In all cases the dipo¬ 
lar symmetry axis is oriented in the rotational equatorial 
plane and the right panels show the resulting pulse pro- 
hles (in arbitrary units) which an observer, also located in 
the rotational equatorial plane, would detect. In all cases 
the core temperature is the same but the star’s distance 
has been adjusted to give the same average hux (see Fig. 
8). Number pairs within parentheses on the left panels 
give (TavejTeff), as dehned in Eq. 15 and 16, resp., in units 
of 10® K. All 5 models have almost the same maximum 
surface temperature Tmax — 8.45 x 10® K but different 
minimal temperatures Tmin- 

Fig. 8: Observable spectra for the hve surface tempera¬ 
ture distributions and observable pulse prohles, labelled as 
“a” to “e”, shown in Fig. 7. The stars, with radii R = 11.4 
km and Roo = 14.28 km for a 1.4 Mq, are assumed to be 
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at distances of 100, 142, 131, 202, and 220 pc, resp., to 
produce almost identical observable spectra in the X-ray 
band {Njj = 1 x 10^° cm“^ for interstellar absorption) 
but resulting in significantly different fluxes in the optical 
range. 

Fig. 9 : Distribution of the surface temperature in terms of 
the percentage of the superficial area whose temperature 
is lower than T, for the five models of Fig. 7 and 8 
Fig. 10: Fit of the spectrum of RX J1856.5-3754. Dotted 
lines show the two blackbodies fit to the data from 
Triimper et al. II2()()4II . The continuous line show our re¬ 
sults: the star has a radius R = 14.4 km and R^o = 17.06 
km for a 1.4 M©, at a distance of 122 pcs [Nh = 1.6 x 10^° 
cm“^ for interstellar absorption) and the observer is as¬ 
sumed to be aligned with the rotation axis. The magnetic 
field structure corresponds to model c of Figure 6 adjusted 
to the 14.4 km radius with T}, = 6.8 x 10^ K, resulting in 
= 4.62 X 10*5 K and = 8.54 x lO'^ K. 
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